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Abstract 

Nonlinear dynamical systems possessing an invariant subspace can display in- 
teresting dynamical behavior, such as on-off intermittency and bubbling. This 
letter shows that a class of such systems have amazing features of (1) supersen- 
sitivity to small input signals and (2) robustness of the supersensitivity in the 
presence of noise. These features make the systems very promising as small 
signal detectors. 
PACS number(s): 05.45. +b; 



Nonlinear dynamical systems with an invariant subspace due to symmetries or other 
constraints are of great interest. An example is the synchronization of coupled chaotic sys- 
tems 0. Such systems can display interesting and unusual dynamical behaviors, such as 
on-off intermittency and bubbling ||. In on-off intermittency, the invariant manifold is 
slightly unstable, and the system can remain close to the invariant manifold for long periods 
of time, interrupted only by some occasional large bursts away from the invariant manifold. 
In bubbling the invariant manifold is stable. However, there are unstable invariant sets em- 
bedded in the chaotic sets, and small perturbations of noise or parameter mismatches can 
result in large intermittent bursts. This can be harmful in connection with applications of 
synchronization, such as in secure communication Q, because high-quality synchronization 
is destroyed by bubbling ||. 

The purpose of this letter is to show that such systems are very sensitive to small constant 
or time-dependent input signal. With an additional symmetrical condition, the sensitivity is 
robust to external noise, which makes the systems very promising for potential application 
in weak signal detection. 

Inspite of the variety of such systems, their behavior often can be described by the fol- 
lowing equations: 

y n+ i = G(x n ,y n ,a), (1) 
x n+1 = F(x n ,y n ,a), (2) 

where G(x n , 0, a) = 0, and the variables y n and x n represent the distance from the invariant 
manifold y = and the dynamics within the invariant subspace, respectively. In general, 
y and x are vectors. Note that many generic properties of the above phenomena can be 
observed in some very simple systems, and we consider x and y as one-dimensional variables. 
Here a is a parameter which may change the dynamics within the invariant subspace as well 
as their stability. We are interested in the behavior of the system close to the invariant 
manifold where the systems can be represented by the approximate linear dynamics: 

Vn+i = g(xn,a)y n + 0(yl) } (3) 
x n+1 = f(x n ,a) + 0(y n ). (4) 

The nonlinearity of the systems serves to keep the solution bounded. Usually, chaotic sig- 
nals have quickly (exponentially) decaying correlation. It is often plausible to assume that 
the chaotic signals are uncorrelated when considering the long time behavior of the system 
in Eq. (0). Based on this assumption, and for the sake of simplicity and without loss of 
generality, we are led to consider the behavior of the following simple random driven map 

y n+1 = ax n y n + 0{y 2 n ), (5) 

where x n is a random driving signal. Introducing the variable z n = In \y n \, Eq. fl5|) becomes 

z n+ i = z„ + ln|x n | +ln|a|. (6) 

The critical value of the parameter a is defined by In |a c | + (ln \ x n \) = 0, where (• • •) represents 
time average. Close to the critical point 5 = (a — a c )/a c <C 1, Eq. (H) can be rewritten as 

z n +i = z n + S + f n , (7) 



where £ n = In \x n \ — (In \x n \) is a random variable with vanishing mean and variance D. 

To analyze the long time behavior, map (|7|) can be replaced by the corresponding stochas- 
tic differential equation, which is the equation for Brownian motion in one dimension with a 
drift. The corresponding Fokker-Planck equation is 

dW dW Dd 2 W 

and the static solution of the probability density W is W(z) = C exp(az) with a = 25/D. 
In the variable y, the distribution becomes W(y) = C\y\ a ~ l . 

Now let us consider that there is a small positive constant input p to the system of Eq. (^), 
namely, y n+ i = ax n y n + 0(y 2 ) +p, where p is the order of lCT m , m ^> 1 (suppose the maximal 
value of y n has the order of unity). For simplicity, we suppose that y n > for initial value 
y > and p > 0. For y n ^> p, the behavior of the system is governed by Eqs. (Q) and (0). 
The effect of the small input can be regarded as a reflecting barrier to the Brownian motion 
of the system ([?]), i.e. z > —m. On the other hand, the state of the system is bounded by 
the nonlinearity of the system. We can introduce a parameter r to represent the effect of 
up-boundary of the system. Based on these considerations, the behavior of the system with 
input p can be understood by the Brownian motion confined between the two boundaries. 
The property of the system is determined by the competition between the constant drift 
5 and the diffusion D. The diffusion is dominant for 5 ~ 0, i.e. for parameter a close to 
the critical point a c where the system displays on-off intermittency (bubbling) and can both 
access to the lower and upper boundaries frequently, becoming sensitive to the small input 
and producing large bursts. Otherwise, the drift becomes dominant when a is far away from 
the critical point; for a < a c , the system will spend most of time close to the lower boundary 
and produce rare large bursts, and the small input does not lead to significant large output in 
the system; for a > a c , the system will spend most of time close to the upper boundary and 
access to the level of the input rarely, and the small input does not have significant effects 
on the system behavior also. We can expect that the system is sensitive to small input when 
it is on-off intermittent. 

The above consideration leads to the normalization condition JI Cy a ~ 1 dy = 1, which 
gives C = a/(r a — p a ). Now we can estimate the amplitude of the output signals by the 
ensemble average 

(y) = C7 f><fr= I T ^" P , (9) 
Jp I + a p — 1 

where (3 = {t/p) 01 . 

If (3 ~ 1, the small input can change the behavior of the system greatly. For the conditions 
|a| <S 1, r 3> p and \a\ ln(r/p) 1, one has j3 ~ 1 + a ln(r/p), and 

(y) « rr^Tv ( 10 ) 

ln(r/p) 

Eq. (|10D shows that the average value decreases to zero with the decrease of input p only 
logarithmically, suggesting that close to the critical point, a very small input p can produce 
a relatively large output, i.e., the system is supersensitive to small input. A measure of the 
sensitivity can be 



For example, with r = 1 and p = 10~ 15 , the value of S is about 2.9 x 10 13 . 

To demonstrate the validity of the above analysis, we carry out simulations with the 
following two systems of the form y n+ \ = ax n f(y) +p. For system I, f(y) is a piecewise linear 
map 

f %{-c\ -c 2 -y), y < -ci, 
f(y) = y, \y\ < (12) 

[ ^(c x + c 2 -y), y > ci, 

where the parameters ci and c 2 are chosen so that y n > for the positive initial value yo and 
p, i.e., the bursting in the system is symmetry breaking^. We use ci = 1 and c 2 = 2 in our 
simulations, with x n uniform on [0, 1] and thus a c = e = 2.71828 • • • and D = 1. For system II, 
f(y) = sin(y), and x n is a chaotic signal generated by the logistic map x n+ i = 3.75x n (l — x n ) 
which gives x n a distribution with singularities, and a c « 1.673 and D = 0.2. S is estimated 
closed to the critical point for these two systems. With p = 10~ m , S as a function of m is 
shown in Fig. 1(a). The analytical estimation in Eq. ( |Tl"{ ) with r = 1.8 and r = 1.4 gives good 
approximation to the simulation results. Fig. 1(b) shows the dependence of S on parameter 
deviation 5 from the critical point. The sensitivity is maintained over a large range of the 
parameter a. 

The feature of supersensitivity is maintained even for time-dependent signals, as for ex- 
ample 

/ p, 0<n<N/2, 

Pn + N = Pn={ _ pj N/2 <n<N. (13) 

In order for the system to have symmetrical response property to positive and negative inputs, 
we require that the map f(y) have odd symmetry f(—y) = —f(y), and on-off intermittency 
in the system is symmetry-breaking, so that a positive (negative) small input will eventually 
lead to only positive (negative) output. For iV 3> iVo, where No is the relaxation time of the 
system, the sensitivity can be measured by Eq. flTT|). 

The feature of supersensitivity makes the systems very promising for application as sen- 
sitive device for small signals. In the context of application, we should consider the behavior 
of the system in the presence of additive noise, namely, y n+ \ = ax n f(y n ) + p n + e n , where 
e n is a small Gaussian white noise with zero mean and standard deviation a . We can study 
the long time behavior of the system by the corresponding stochastic differential equation 
dy/dt = (5 + £)y + p + e. The Fokker-Planck equation is 

f = -|{[(^f)»+^} + iJ[W + on a*, 

and the static solution under the adiabatic condition N ^ N is given by 



W(y) = C(y 2 + ^)( a - 1 )/ 2 exp[^arctan^]. (15) 



a 



This distribution, however, is very complicated for evaluating (y). To simplify the calculation, 
we employ the similar heuristic boundary conditions in the above. Under the conditions 
a v^D, p ~ cr, we approximate the distribution by 

w(y) = I ° lyrl ex ^vk s ^ n ^ \y\ * *>> m 

y 1 \ 0, \y\ < p. y 1 



It is clear that by the limit a — > 0, we come back to the result for the noise-free case. With 
this approximation, we obtain that close to the critical point, 



(y) 



j: T yW(y)dy 




(17) 



r_ T W{y)dy ln(r/p) 



where R = p/a provides a natural measure of the signal-to-noise ratio. 

The above analysis is demonstrated by numerical simulations in the presence of noise. 
Fig. 2(a) is a typical response of the system to a noisy small signal and Fig. 2(b) shows (y n ) 
for different values of R. Fig. 2(c) displays the dependence of (y) on R. It is seen that the 
above approximate analysis gives a good account for the results in a large range of R. Over 
a wide range of R, (y) is very close to that of the noise-free case. The supersensitivity is 
thus robust to additive noise. This feature of sensitivity is quite different from that of the 
sensitivity near the onset of a period-doubling bifurcation in many dynamical systems ||. 
There the system is only sensitive to perturbations near half the fundamental frequency of 
the system for bifurcation parameter very close to the onset point. 

To examine the performance of the system as a small signal detector, we calculate the 
probability of bit error P& in the presence of additive noise. The detection is done by 
examining the time average of the output y n in the duration of a input bit bk, namely 



variable Sk is expected to fluctuate around (y). Although y n cannot assumed to be uncorre- 
cted, for very large N, it might still be plausible to assume that Sk approaches a Gaussian 
distribution with an average (y) and a variance Dn = A/N, especially in the case that R 
is small and y n has comparable distribution to positive and negative values. Based on this 
assumption P b can be evaluated approximately as 



where (y) is given by Eq. flTfl). 

In our simulations, we estimate with 10 6 random bits in the system I at a = 2.6 for 
input levels p = 10~ 4 and p = 10 -6 . The quantities N Q and A are estimated in simulation 
with constant input, giving iVo = 350, A = 5 for p = 10 -4 and N = 700, A = 4 for p = 10~ 6 . 
Both the results of Pb from simulations and from estimation in Eq. (|i~8[) are shown in Fig. 3 
for N = 3iV , 5N , 10iV , 15iV . The parameter r used to fit Eq. ([Tj|) to the simulation results 
is r = 1.4 for p = 10~ 4 and r = 1.1 for p = 10~ 6 . It is seen that the estimation can be 
quite good for N much larger than Nq. For iV comparable to iVo, the effects of the transient 
process during the relaxation time becomes significant, and the estimation deviates from the 
simulation results. For R getting larger, the bursting behavior become more asymmetrical, 
and Sk can no longer be approximated by a Gaussian distribution and the estimation also 
deviates from the simulations results. 

The above results show that detection error can be quite low even for a small signal with a 
level much lower than the environment noise if the signal has a bit duration much larger than 
the relaxation time of the system. Related to this, there seems to be a frequency cutoff above 
which detection becomes unreliable. This frequency gets larger as the input level increases, 
because the relaxation time becomes shorter for higher level of input. 

The above properties of supersensitivity and its robustness in the presence of noise is 
universal in a general class of coupled symmetrical systems displaying on-off intermittency 





(18) 



with symmetry-breaking. The sensitivity is due to the power-law distribution of on-off inter- 
mittency of y in a wide interval 10 _m < \y\ < r, so that the system can both access to the 
level of the small input and at the same time produce frequent large bursts. If the maps are 
not coupled to random or chaotic driving x n , there is no diffusion in the system (D = 0), 
and the small input does not have significant effect on the system output. For uncoupled 
non-on-off maps, namely, y n +\ = afiUn), if a < 1 (for the maps in the above), the fixed 
point y = is stable, and the small input cannot produce large output at all; if a > 1, the 
state y n can no longer come to the level of a small input of the order p = 10~ m (m > 1) 
with significant frequency, and the output will not manifest the small input. In both cases, 
the systems do not possess the sensitivity in the coupled, on-off intermittent systems. The 
symmetry-breaking of the bursting also plays an important role in the sensitivity and ro- 
bustness, because under this condition, a transition of the state between y > and y < is 
determined only by the switch of the small signal. If the bursting is not symmetry-breaking, 
there are additional transitions between y > and y < induced by bursting states, which 
will degrade the sensitivity and robustness. 

In conclusion, we demonstrate that a class of nonlinear dynamical systems having an 
invariant subspace and displaying on-off intermittency and bubbling have s the feature of 
supersensitivity to small constant or time-dependent input signals. With an additional odd 
symmetry condition, the sensitivity is robust to additive noise. The features make the systems 
very promising for useful application as sensitive devices. 
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Figure Captions 



Fig. 1. (a) Dependence of the sensitivity S close to the critical point of the systems on 
input p = lCT m . (b) Dependence of S on 6 for p = 1CT 10 . 

Fig. 2. (a) An example of the bursting behavior of the system I with a = 2.6, N = 4000, 
p = 10~ 5 and R — 0.1. (b)Time series of the ensemble average (y n ) over 5000 samples 
for the system I with a = 2.6, p = 10~ 5 and iV = 4000. The three plots are: (1) for 
noise-free case, (2) for R = 0.2, and (3) for R = 0.05. (c) Ensemble average (y) close 
to the critical point as a function of R for constant input p = 10~ 5 . The solid lines are 
estimation of Eq. (|17|). 

Fig. 3. The probability of bit error as a function of R for different levels of input 
and different bit durations. The solid lines are estimation of Eq. (a) p = 10~ 4 , 

iVo = 350, A = 5 and r = 1.4. (b) p = 10" 6 , iVo = 700, A = 4 and r = 1.1 
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